Topological edge states at a tilt boundary in gated multi-layer graphene 
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Despite much interest in engineering new topological surface (edge) states using structural defects, 
such topological surface states have not been observed yet. We show that recently imaged tilt 
boundaries in gated multi- layer graphene should support topologically protected gap less edge states. 
We approach the problem from two perspectives: the microscopic perspective of a tight-binding 
model and an ab-initio calculation on a bilayer, and the symmetry protected topological (SPT) 
states perspective for a general multi-layer. Hence we establish the tilt boundary edge states as 
the first concrete example of edge states of symmetry enriched Z-type SPT, protected by no valley 
mixing, electron number conservation, and time reversal T symmetries. Further we discuss possible 
phase transitions between distinct SPT's upon symmetry changes. Combined with recently imaged 
tilt boundary network, our findings offer a natural explanation for the long standing puzzle of sub-gap 
conductance in gated bilayer graphene, which can be tested through future transport experiments 
on tilt boundaries. In particular, the tilt boundaries offer an opportunity for in-situ imaging of 
topological edge transport 



I. INTRODUCTION 

Graphene has garnered interest from broad spectrum 
of communities, ranging from those aiming at atomic 
scale circuit devices to those searching for new topo- 
logical phases. Both communities sought after ways to 
gap the massless Dirac spectrum. The realization of 
a gate-induced band-gap in the Bernal stacked bi-layer 
graphene [1] following the prediction in Ref. [2] brought 
the holy grail of graphene based transistor one step closer 
to reality. However, the sub-gap conductance measured 
by Oostinga et al. [1] with weak temperature depen- 
dence well below the optically measured gap as large 
as 250 meV[3] introduced a new puzzle and obstacle: 
the gapped bilayer is not as insulating as it should be. 
Dominant transport along physical edge of the samples 
proposed earlier by Li et al. [4] have been ruled out by 
Corbino geometry measurements [5], which observed two- 
dimensional variable range hopping type temperature de- 
pendence, independent of geometry. In this paper we 
predict existence of topological gapless channel of trans- 
port along recently imaged AB-BA tilt boundary network 
[6-8] which solves the puzzle. 

The predicted topological edge state holds the promise 
of the first realization of topological surface (edge) state 
hosted by structural topological defect. Though there 
has been much theoretical interest in topological gap- 
less modes hosted by structural topological defects [9-11] 
no such topological gapless mode has been observed so 
far. The lattice dislocations in three dimensional crystals 
previously discussed occur deep in the sample that is not 
directly accessible. However, the tilt boundary of inter- 
est have recently been observed [6-8]. The tilt boundary 
is a structural topological line defect along which each 
neighboring layer is displaced by one inter-atomic spac- 
ing. Such defect can occur due to the third dimension 
added by the stacking of the graphene layers; it forms a 
boundary between two inequivalent stacking structures 
frequently referred to as AB and BA. Here we show that 
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FIG. 1. A typical AB-BA tilt boundary under strain. The 
blue (red) filled circles mark the a (b) sublattice sites, jij 
represent hopping matrices for a tight-binding model. 



the tilt boundaries host gapless modes of topological ori- 
gin and form the first example of a naked structural de- 
fect hosting topological electronic states. 

Topological aspects of gapped multi-layer graphene 
have been previously discussed [12] and it was pointed 
out that they should exhibit quantum valley Hall effect 
with corresponding edge states. However, to this date 
there has been no experimental detection of proposed 
edge state [4, 13]. Moreover, little is known about how 
the topological aspects of gapped multi-layer graphene 
relates to topological insulators [14, 15]. The idea of 
classifying different topological insulator (superconduc- 
tor) candidates based on symmetries [16] have played a 
key role in the field of topological insulators. In par- 
ticular the observation that additional symmetries such 
as the crystalline symmetries can enlarge the possibili- 
ties of topological phases [17, 18] led to the discovery of 
three-dimensional topological crystalline insulators [19]. 
On one hand we propose feasible experiments to detect 
topological edge states at naturally occurring tilt bound- 
aries. At the same time, we make first concrete applica- 
tion of the SPT approach [17] for two dimensional (2D) 



system and study a large class of gapped graphene sys- 
tems placing the quantum valley Hall insulator in the 
larger context and predicting conditions for topological 
superconductors. 

The rest of the paper is organized as follows. In sec- 
tion II we show that a AB-BA tilt boundary in gated 
bilayer graphene supports gapless edge states through ex- 
plicit microscopic calculations. Specifically we consider 
an abrupt boundary in tight-binding model and then in- 
vestigate the effect of strain using ab-initio calculation. 
In section III we show that these edge-states are pro- 
tected by no valley mixing, electron number conservation, 
and time reversal (T) symmetries within the framework 
of SPT. Hence we identify chirally stacked gated A/'-layer 
graphene layers as time-reversal symmetric Z-type SPT. 
In section IV we discuss experimental implications. Fi- 
nally in section V we summarize the results and comment 
on practical implications. 

II. BI-LAYER TILT BOUNDARY EDGE STATE 

Fig. 1 and Fig. 2(a) show tilt boundaries of interest 
in gapped Bernal stacked bi-layer graphene. In the case 
sketched, strain is concentrated at the tilt boundary with 
the top layer stretched by one inter-atomic spacing with 
respect to the bottom layer. For a general orientation, 
tilt boundaries can involve both strain and shear. As 
the tilt boundaries in layered graphene form a type of 
topological line defects in structure, they can be charac- 
terized using the tangent vector t and the Burger's vector 
b. The tangent vector t points along the tilt boundary 
which can point along any direction with respect to the 
Burger's vector b. When the tilt boundary only involves 
strain as in the case depicted in Fig. 1 and Fig. 2(a), the 
b is perpendicular to t. In the opposite extreme limit of 
6 II t, shear is concentrated at the boundary. Indepen- 
dent of the angle between b and t, the Burger's vector 
magnitude is the inter- atomic spacing i.e. \b\ = a for 
a bilayer system, as it is shown explicitly for the strain 
tilt boundary in Fig. 2(a). Since \b\ is a fraction of the 
Bravis lattice primitive vector magnitude >/3a, the bi- 
layer domain boundary is a partial dislocation from quasi 
two-dimensional view. In a general mult-layer a vertical 
array of these partial dislocations form a tilt-boundary. 
In typical samples, the domain wall separating the AB 
and BA stacked domains has substantial width spanning 
5-20 nm and the angle between b and t ranges between 
0° and 90° [6]. 

Fig. 1 allow us to makes two important microscopic 
observations about tilt-boundaries. (1) As the bound- 
ary requires a shift of one-layer with respect to the other 
by one inter-atomic spacing along the bond direction, 
there are three natural directions for the tilt-boundary 
to run for each fixed angle between b and t. (2) The 
boundary is arm-chair for b \\ t (pure shear) whereas it 
is zigzag for b -L t (pure strain case shown in Fig. 1). 



2 




-0.5 0.5 

k {2n/c) 
(b) 



AB BA 

(c) 

FIG. 2. A tight binding study in the presence of abrupt AB- 
BA tilt boundaries in gated bilayer graphene. (a) Schematic 
representation of the domain wall under strain. The grey 
(black) lines denote the upper (lower) layer, and blue (red) 
solid circles denote the A (B) sublattice points. The tilt 
boundary along tangent vector t \\ y is shaded. The Burger's 
vector b = ax (green arrow) for interatomic spacing a accounts 
for the difference between the green solid line and the green 
dashed line, (b) The resulting band structure. Edge states 
are marked in magenta, (c) Schematics of valley- momentum 
locked edge states at a tilt boundary. 



Based on observation (1) we expect a given type of tilt 
boundary to possibly form a triangular network seen in 
experiments [6-8]. The observation (2) combined with 
earlier microscopic studies of boundary condition effects 
on edge states in Ref. [20] implies that electronic spec- 
trum at tilt-boundaries with b \\ t will be gapped though 
the gap magnitude will be small when the tilt boundary 
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is spread over finite width. 

In this section, we consider the electronic structure of 
tilt boundaries with b ±t and return to more general case 
in the section IV. As it is shown schematically in Fig. 1, 
the bulk of each domain is gapped in the presence of inter- 
layer hopping and the external electric field. The latter is 
important for breaking the inversion symmetry between 
the layers and gapping otherwise touching bands [1, 2]. 
Below we present two separate microscopic calculations 
of the AB-BA tilt boundary electronic structure for 6 ± 
which shows gapless edge states. 



A. Tight-binding model 

We consider a tight-binding Hamiltonian with nearest- 
neighbor intra- and inter-layer hopping. For the AB- 
stacking region (see Fig. 2(a)), 

2 

Hab =-t J2 E (^™," + + ^™!n-l) 

-t±j2^ll]Jb^^l + h.c., (1) 

m,n 

where i = 1, 2 is a layer index, t and t± are intra-layer and 
interlayer nearest neighbor hopping respectively, and A 
is the chemical potential difference between two layers 
due to the gate voltage. (m,n) labels the position of 
the two site unit-cell with am^n and b^^n annihilating 
electrons at the two sites of layer i. For the BA stacked 
region, the only change is in the inter-layer term with 

{—t_\_ bm]n aln]n + h.c) replacing the inter-layer term 

m,n 

in HaB' As we address the effect of strain through ab- 
initio simulation, we focus here on a sharp tilt boundary 
as shown in Fig. 2. 

We plot the energy spectrum in Fig. 2(b) with the 
model parameters set to be t = 2.8eV, t± = 0.4eV, and 
A = 0.5eV. The size of the system was 200 unit-cell in 
each direction under periodic boundary condition with 
each domain spanning 100 unit-cell width separated by 
two sharp tilt boundaries. From the spectra it is clear 
that K and valleys each have two edge states per spin. 
Further investigation of the wave function shows right 
and left moving edge states associated with the given 
valley are spatially separated between the two edges: the 
edges offer valley filtering (see Fig. 2(c)). 



B. Ab Initio Simulation 

The electronic structure in realistic tilt boundary will 
be affected by both the span over which the lattice struc- 
ture transition from AB to BA stacking, and the strain 



concentrated at the tilt boundary. To address these is- 
sues we further carried out first-principles calculation 
using density functional theory (DFT) within the local 
density approximation (LDA) [21, 22]. We constructed 
a periodic supercell with two tilt boundaries character- 
ized by Burgers vectors identical to the ones considered 
within tight binding calculation, but the domain wall is 
set to have finite width over which the one unit-cell mis- 
match is spread. We have tried several configurations 
with domains and tilt boundaries of different widths to 
find qualitatively similar results. In the rest of this paper 
we focus on a representative example with 3.1 nm wide 
domain wall between 1.5 nm wide domains. The choice 
of narrow width for the domains was due to the limita- 
tion in the simulation capability. To this 2-D system, we 
then applied a slightly exaggerative perpendicular elec- 
tric field of 5 V/nm. By relaxing the domain boundaries 
until the forces reaches below 0.5 eV/nm, we took the ef- 
fect of both strain and the width into account. We have 
tried several configurations with domains and tilt bound- 
aries of different width but otherwise similar setting to 
find qualitatively similar results. 

The DFT-LDA simulation results are presented in Fig. 
3, which confirms the existence of the gapless edge states 
predicted by the tight-binding model. The electronic 
structure along the extended direction of edge is shown 
in Fig. 3(a), focusing on the region near K point. In this 
figure two gapless ID Dirac dispersion is clearly resolved 
from the gapped bulk specta, with two distinct Dirac 
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FIG. 3. (a) The energy bands of AB-BA bilayer graphene 
connected two domain walls depicted in Fig. 2. The Fermi 
level is indicated with the blue dashed line and the bands 
with marked linear dispersion relation intersecting near the 
Fermi level are outlined by the purple (light) curves, (b) Top 
and side views of the charge distribution for a state near the 
Ki point of (a). The yellow rectangles indicate the AB-BA 
domain wall structure. 
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points Ki and K2 in the vicinity of the K point. The 
strain concentrated at the tilt boundary causes energy 
sphtting of Ki and K2 states; the energy of ID Dirac 
points is increased by compressive strain and decreased 
by tensile strain. This energy splitting will become neg- 
ligibly small in realistic tilt boundaries with much wider 
span, as the strain will become smaller. Fig. 3(a) shows 
that the existence of gapless edge states found in our 
tight-binding calculation are robust against long range 
perturbations such as tilt boundary width or subtle bond- 
length variation inside the domain wall, as well as the 
strain at the tilt boundary. 

We now turn to the spatial distribution of the gapless 
edge states. For illustration. Fig. 3 (b) shows the wave- 
function amplitude of the state slightly above the Ki 
point. As expected from the tight-binding results, the 
charge distribution is prominent inside the domain wall 
but rapidly decays away from the tilt boundary. The 
Dirac points of edge states and the decaying feature of 
edge states always exist, making us believe it is a uni- 
versal feature. Meanwhile, the charge for the states near 
Ki is highly localized on the layer subject to compressive 
strain, regardless of the direction of the applied electric 
field. The situation for the states in vicinity of K2 is simi- 
lar except that the charge prefer to highly localized on the 
layer under tensile strain. Hence forming a layer selective 
contact to an isolated edge state could be a mechanism 
for valley filtering[12, 23]. 

Finally, we comment on the so far ignored effect of in- 
teraction. If the edge boundaries have substantial width, 
forward scattering part of the Coulomb interaction will 
be the dominant correlation effect to the edge states and 
lead to Luttinger liquid behavior [24]. 



III. TILT BOUNDARY EDGE STATES AND 
Z-SPT 

In order to address the robustness of the edge states, 
we investigate topological aspects of the low energy effec- 
tive theory in the continuum limit. We first show that the 
AB-BA tilt boundary can be mapped to a gate-polarity 
boundary of uniform bi-layer. Based on this mapping and 
results of Refs. [13, 20] on the gate-polarity boundary, we 
discuss the valley Chern number of the tilt boundary edge 
states. We than apply the notion of SPT [17] and iden- 
tify chirally-stacked multi-layer graphene as a realization 
of Z-type SPT, protected by time- reversal (T), absence 
of valley mixing, charge conservation symmetries. This 
identification enables us to address effects of symmetry 
changes: topological quantum phase transitions. There 
are recent studies of such perturbations for specific cases 
such as Rashba spin-orbit coupling [25] and magnetic or- 
dering [26] . Through our first application of SPT classifi- 
cation scheme by Wen [17] to a concrete physical system 
of multi-layer graphene, we obtain an exhaustive system- 
atic study of topological quantum phase transition pos- 
sibilities. 



A. Valley Chern Number 

The low energy effective Hamiltonian near the K valley 
for uniformly AB or BA stacked bi-layer is 

AB / BA 

Hj^ =VFkx{fix^Xo)^VFky{fiy^Xo) (2) 
-(A/2)(/io A^) + y (/i^ A^) =F y (My ^ A^), 

where — or + sign should be used for AB or BA stacking 
respectively. In Eq. (2), vp = 3ta/2, /i^'s are Pauli ma- 
trices acting on the sub-lattice indices, and A^'s are Pauli 
matrices acting on the layer indices. The effective Hamil- 
tonian near K' is Hk' (kx^ky) = Hk {—kx^ky). Now it is 
straight forward to show that BA stacking is equivalent 
to AB stacking subject to the opposite gate polarity. At 
zero field the Hamiltonian of the AB stacked bilayer can 
be transformed to that of the BA stacking by interchang- 
ing the two layers via the following unitary transforma- 
tion: H S'^HS^ with S = fiQ^Xx- For a gated bi-layer 
however, the gate polarity has to flip since 

Ho <S> X,H^^{A)txo ® A, = H§^{-A). (3) 

Hence at the level of low energy effective theory, the tilt 
boundary between AB and BA stacking under uniform 
external field is equivalent to the gate polarity domain 
wall of structurally uniform bi-layer proposed by Martin 
et al [13]. 

The above equivalence combined with earlier results on 
valley Chern number of gated chirally stacked multi-layer 
offers the topological origin of the helical edge states ob- 
served in the microscopic calculation of section II. First 
for bi-layer and then for general TV-layers, it was shown 
that low energy effective theory of chiral stacked A^-layer 
under uniform vertical electric field A has finite Chern 
number per spin for each valley of equal magnitude and 
opposite sign [13, 20, 27, 28]: 

TV 

Ck = -C'^ = ySgn(t^A). (4) 

The Chern numbers in Eq. (4) can be obtained by inte- 
grating the Berry curvature over momenta {kx^ky) con- 
tinuing the linearized dispersion to infinity. This com- 
bined with the equivalence relation of Eq. (3) means the 
valley Chern numbers change sign at the tilt boundary. 
Such sign change leads to AC^ = —ACr = Nsgn{t±A) 
across the tilt boundary and N branches of valley heli- 
cal edge states [20], as long as the two valleys K and 
remain distinct. Hence the two valley helical edge states 
per spin observed in section II originate from the val- 
ley Chern number change across the tilt boundary, as in 
the gate polarity boundary edge states [13, 20]. Hence, 
our prediction is the tilt boundaries will be the first ex- 
perimentally observed crystalline topological defects to 
host topological gapless mode due to change in the Chern 
number. 
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B. Chiral Multi-layer as Z-type SPT 

We now apply the procedure for identifying the class 
of SPT based on symmetries of free fermion Hamiltonian 
developed by Wen [17], which predicts possible number 
of protected edge (surface) states. This procedure allows 
us to consider additional symmetries in the multi-layer 
graphene in addition to the C, T, and P taken into ac- 
count in the pioneering work by Schnyder et al. [16], and 
by Kitaev [29]. The procedure consists of three steps: 

(1) Find a gapless Dirac Hamiltonian (by keeping the 
kinetic term only) with the same symmetries. Then we 
find all the symmetry preserving mass terms that can 
gap out the gapless part and are amenable to classifi- 
cation using Clifford algebra. This is based on the as- 
sumption that the SPT order is robust as long as the 
energy gap stays finite and the symmetries remain the 
same and hence any gapped Hamiltonian can be adia- 
batically transformed into a gapped Dirac Hamiltonian. 

(2) Express the Hamiltonian and the conserved quan- 
tities associated with symmetries in the Majorana basis. 
This leads to the Clifford algebra (i.e. real representation 
of the Dirac algebra) associated with the gapless part of 
the Hamiltonian. (3) Find the space of mass matrices 
that ant i- commutes with all the generators of this Clif- 
ford algebra. The resulting space may have disconnected 
pieces, the number of which gives the classification of 
the SPT. Two mass matrices are topologically distinct if 
and only if they belong to two different pieces. Applying 
this procedure to chiral multi-layer graphene will enable 
us to study phase transitions into different SPTs upon 
symmetry changes. 

For chiral multi-layer graphene system, we assume no 
valley mixing, electron number conservation ([/(l)c), and 
time reversal (T) symmetries. The relevant gapless Dirac 
Hamiltonian is: 



H = wy j d?x (^) ^p^d^ + p^dy) ^ {x) , 



(5) 



where = {ca^k.ca^k' ,cb^k,cb^k'), and pi = iix ® 
Tz ^ Inxn, and p2 = (8) To (8) /nxn, in which n is given 
by the number of layers. This Hamiltonian can be writ- 
ten in the Majorana fermion basis using the following 
decomposition 



t-/Lt,r (•^) 



7+,M,T (x) + il-,ix,T (x) 



(6) 



where 7 denotes the Majorana fermion satisfying 

1 , {7+,^,T,7-,M,r} = 0, (7) 



7^ 



where fi (r) denotes the A, or B sublattice {K or 
valley) indices, and a denotes the flavor of the Majorana 
fermions (+ or -). In the Majorana fermion basis the 
Hamiltonian is represented as follows: 



H 



■'I 



(fx 77 (x) Aq {x) , 



(8) 



where ^ is a real anti-symmetric matrix (differential op- 
erator), and rj is an eight component vector whose com- 
ponents are 'ya,fi,T- 

Now we express the conserved quantities associated 
with the symmetries of the Hamiltonian in the Majorana 
fermion basis. First, no valley mixing combined with to- 
tal electron number conservation symmetry leads to sep- 
arate conservation of the electron number at each valley 
Nx and Nx'- Hence the total electron number Nc = 
Nk + Nk' , and the valley polarization Ny = Nk — Nk' 
are conserved. In the Majorana fermion basis, 

Ncy = ^ J V (^) QcyV (^) 7 with 

Qc = iOLy ^ fio^To^ I, Qv = i<^y ^ Po^Tz ^ I, (9) 

where ai (r^) Pauli matrices act on the Majorana flavors 
(valley indices). So deflned Qc and Qv satisfy Q^ = 

Ql = -1. 

Under time reversal symmetry T, K and valley in- 
dices are exchanged i.e., T : c^^k ^ c^,k' - Hence T acts 
like Tx in the valley basis with the matrix part of the time 
reversal operator satisfying = 1. On the other hand, 
no valley mixing implies the Hamiltonian is invariant un- 
der (c^,k,c^,kO {c^x,Ki —Cfx,K') transformation which 
acts like in the valley basis. Hence, in the presence of 
no valley mixing symmetry, we can define a new time re- 
versal operator 6 : (c^,k,c^,kO ^ {c^,k' , -c^,k), which 
acts like r^T = iry. In terms of Majorana fermions 



Q = fjio^ {iTy) (8) / , 



(10) 



In order to find the relevant Clifford algebra, we need 
to form anti-commuting generators in terms of A in Eq. 
(8) combined with symmetries, Qc, Qv, and Q. How- 



ever, the symmetries require 



w4., Qc 



A, Qi 







and 1^,6 1=0. Moreover, symmetry operators satisfy 

Qc& QQc, Qy© and QcQv QvQc rela- 

tions. Using these relations, we find the full set of gener- 
ators of the relevant Clifford algebra as pi, p2, and 



P3 = QQvQc 



P4 = e , P5 = ©Qy, (11) 



as {A, Pi} = for i = 3,4,5. The resulting full set of 
anti-commutation relations is 

{pi, Pj} = "^Qij , 9i,j = diag (1, 1, 1, -1, -1) , (12) 

and it defines a Clifford algebra Cliff(3,2). 

Now, we will find the space of mass matrices. Cm, 
that can gap out Dirac Hamiltonian associated with 
this Clifford algebra in order to obtain SPT classifica- 
tion. The mass term with matrix representation Hm = 
i Si J Mi,jViVJ should satisfy the following algebra: 



Mpi = -piM, 



-1, 



(13) 



where we have normalized M. Solving the above equa- 
tion yields the allowed space for the mass matrix. Cm- It 
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has been shown [17] that Cm which solves Eq. (14) for 
the case of Chff(3, 2) is 

= 1™ U ^ ^n} (^4) 

n^oo 0{m) X 0{n -m) 

The SPT classification is then given by the number of dis- 
connected pieces in the space of mass matrix Cm i-e. its 
zeroth homotopy group: ttq (Cm)- Using Eq. (14), it can 
be verified that ttq (Cm) = ^ [17]. Consequently, each 
class of the time reversal invariant multilayer graphene 
in the absence of intervalley scattering is indexed by a 
Z- valued number: in this case the valley Chern number. 

Now we are in a good position to consider symmetry 
changes. Important to note here that spontaneously or- 
dered phases can be considered alongside systems under 
external field, as once a system is deep inside the ordered 
phase it can be treated within mean-field theory. 

We first consider the symmetry reduction possibili- 
ties while maintaining spin degeneracy (see Table I). If 
we only break the time reversal symmetry, the system 
is characterized by two independent topological indices 
{Ck-,Ck')^ hence the classification is given by Z Z. 
We refer to these SPT phases as intra-valley quantum 
(anomalous) Hall (QAH) states. Such phases may be re- 
alized by placing trigonally-strained graphene [30, 31] un- 
der an external magnetic field, as the sign of the pseudo- 
magnetic field caused by strain is opposite for the two 
valleys. Further reducing the symmetry by introducing 
inter- valley scattering leads to inter- valley QAH state in- 
dexed by a single integer Z: the total Chern number [32]. 
Breaking electron number conservation turns the above 
insulators into superconductors. Following the procedure 
above, we obtain the same classification for the topologi- 
cal superconductors in 2D, resulting in topological valley 
superconductor (TVSC) and intra- or inter-valley topo- 
logical superconductors. Table I summarizes all symme- 
try reduction possibilities and their classifications start- 
ing form gated multi-layer graphene. 

Now we consider extending our classification to take 
the electron spin into account as a dynamical degree of 
freedom. This will allow us to consider interaction effects 
at the level of spin ordering. With both spin and valley 
degrees of freedom, there are three symmetry operators 
related to time reversal: Oi oc r^^ exchanges two valleys, 
©2 oc i(jy which acts on the spin indices fiips spin, and 
T = 9iB2 does both. If both 62 and T are broken, 
the classification reduces to that of spinless electrons (see 
upper part of Table II). However, taking any of these 
symmetry operators into account leads to new classes 
(see lower part of Table II) . 

When all ©1, ©2 (and as a result T) and no valley- 
mixing symmetries are imposed the state corresponds to 
the so-called "layer antiferromagnetic" (LAF) phase pre- 
dicted in Refs. [28, 33] and possibly occurring as a ground 
state in bi-layer graphene at neutrality point [34]. In 
this phase the product of spin and valley of edge quasi- 
particles is locked to their momentum. Therefore, one 



may index this state by its spin-valley Chern number 
Csv [28]. In this phase, for each valley, the edge states as- 
sociated with two spins are counter-propagating. As long 
as these two counter-propagating modes do not couple, 
there can be any number of them per each valley, lead- 
ing to 2Cs'ye^/h spin- valley Hall conductivity [12, 28]. 
However there is a form of symmetry allowed coupling 
between a pair of sets of counter propagating edges; this 
coupling can gap out the edge modes [35]. Therefore, 
the number of symmetry protected edge modes for each 
valley is Csv naod 2. Hence, we obtain a Z2 classifica- 
tion for this phase labeled by (—1)^^^ as supposed to Z 
classification which would be implied by Refs. [12, 28]. 

Another interesting possibility is breaking ©i, and ©2, 
while respecting their product T and no-valley mixing 
symmetry. This leads to the quantum spin Hall (QSH) 
phase, in which Ck,\ = —Ck',],^ and Ck',\ = —Ck,], due 
to T symmetry, while there is no constraint on the Cx,t + 
(7^4 . Hence, unlike usual 2D QSH states [14, 15], we 
obtain Z. This is because multiple type of time reversal 
operators can be defined [36] in the presence of no valley 
mixing symmetry. This and other SPT possibilities are 
summarized in Table II. 

With the classification at hand, we return to what it 
implies for the fate of the gapless edge states at the tilt 
boundary of chirally stacked and gated N-layer graphene. 
In general, a topological phase transition between two 
phases within the same class requires the bulk gap to 
close and reopen (as in inter-plateaux transitions in in- 
teger quantum Hall effect). Gapless edge states are gu- 
ranteed at a physical boundary between such two phases. 
The edge states at the AB-BA tilt boundary we estab- 
lished in the section II are examples of such edge states. 
Hence the edge states will remain gapless as long as the 
time (valley) reversal and charge conservation symmetry 
are maintained. 

On the other hand, symmetry change can either yield 
a trivial phase which does not support an edge state or a 
different type of SPT with different type of edge states. 
According to Tables I and II, ruining the no valley mix- 
ing symmetry is the only way to render the system triv- 
ial and gap the edge states. However this requires large 
momentum transfer which generally requires fine tuning 
unless the unit-cell becomes enlarged either for the en- 
tire system [37, 38] or for the edge through arm-chair 
edge, or short-range disorder such as a vacancy breaks 
A-B sub-lattice symmetry [39]. As both rarely occur, we 
anticipate gapless edge states at most tilt boundaries. In 
particular, the natural zigzag boundary formation for the 
tilt boundaries make such edge states more robust than 
the edge states in gate polarity boundaries [13]. 

Among symmetry change possibilities leading to an- 
other SPT, transition from QVH with spin degeneracy 
to LAF phase where spin is a dynamic degree of free- 
dom is of particular interest as LAF is suspected to 
be the ground state of bi-layer graphene near neutral- 
ity point [34]. Upon this phase transition the nature of 
edge states change from spin degenerate valley helical 
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TABLE n. Classification of spinful SPT insulators on multi-layer graphene. Three kinds of time reversal operators are considered 
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on the spin indices fiips spin, and T = B1B2 does both. ^ The upper three rows show that the classification reduces to that of 
spinless fermion if B2 symmetry is absent. The lower rows show new possibilities that emerge upon taking spins into account. 

^ The classification for the superconductors are identical to that for the insulators. 



states (QVH) to spin- valley Hall edge states (LAP). In 
such transitions, the bulk gap has to close and reopen; 
this is indeed seen in the experiment of Velasco et al [34] . 

IV. CONNECTION TO EXPERIMENTS 

In this section, we discuss the experimental implica- 
tions of our findings. Specifically we propose the trans- 
port through the network of tilt boundaries as a solution 
to the long standing mystery of sub-gap transport [1, 5]. 
Purther we propose feasible experiments to test the pro- 
posal. 

In order to discuss the topological transport through 
the network of tilt boundaries observed in Refs. [6-8] 
we should first discuss the effect of the arbitrary an- 
gle between the Burger's vector b and the tangent vec- 
tor t. Microscopic study of tilt boundaries at vari- 
ous angles between b and t will be presented in the 
future [40]. However, it has been known that a single 
domain Bernal stacked gapped bilayer ribbon should sup- 
port gapless edge state for zig-zag edges, but not for arm- 
chair edges [20, 41]. This is because the arm-chair edge 
enlarges the unit-cell along the direction parallel to the 
boundary and makes the projection of K and valley 
identical. However, Jung et al. [20] showed that the po- 
larity boundary edge states only develop barely visible 
gap which is orders of magnitude less then the bulk gap 
even for a sharp boundary, and the gap decreases quickly 
when the polarity boundary becomes smooth. These ar- 



guments apply to our tilt boundary and the edge state 
will develop a small gap when b \\ t. However, given the 
large width of the observed tilt boundaries we expect that 
all straight tilt boundaries will have nearly gapless edge 
states except those with small angles between b and t. 

When the tilt boundary meanders and changes direc- 
tions, likely there will be portions with small gap seg- 
menting the gapless regions and the transport will occur 
through hopping between the gapless regions. The ob- 
served 2D network of such tilt boundaries would yield 2D 
variable range hopping temperature dependence R{T) cx 
exp(To/T)^/^[42, 43] at low temperatures governed by 
the 2D connectivity and the small characteristic gaps of 
gapped regions. This explains the observed temperature 
T dependence of resistance at low temperatures [1, 5]. 

We propose following experiments to test our proposal. 
(1) Pour terminal transport measurements with two of 
the contacts, say contacts 1 and 3, at two ends of a tilt 
boundary. This would yield highly anisotropic transport 
proving dominant transport along the tilt boundary i.e., 
^1,3 ^ ^2,4- (2) Scanning tunneling spectroscopy mea- 
surements of local density of states. This should mea- 
sure a gapless spectrum at the tilt boundary but exhibit 
a gapped spectrum with the gap magnitude of the opti- 
cal gap away from the tilt boundary. (3) Thermoelectric 
imaging. The mid-gap density of state at tilt boundaries 
would appear in scanning thermopower images. Unpub- 
lished thermopower imaging data by Cho et al [44] in- 
deed show a network with local density of state near fermi 
energy, that is reminiscent of the tilt boundary network. 
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(4) Edge current imaging using scanning SQUID which 
can detect magnetic field generated by edge currents. 



V. CONCLUSION 

We showed that spin-degenerate tilt boundaries of 
gated multi-layer graphene support topological gapless 
edge states protected by three symmetries: time (valley) 
reversal, no- valley- mixing, and electron number conser- 
vation. We demonstrated the existence of gapless edge 
states through a tight-binding model calculation and a 
first principal calculation, where the latter took strain 
effects into account. We then addressed the symmetry 
protection of the edge states and consequences of sym- 
metry changes within the framework of SPT [17]. 

The framework of SPT allowed us to place the 2D topo- 
logical phase supporting the edge states, namely QVH, 
among various topological insulator /superconductor 
phases alongside previously postulated QAH, LAF and 
QSH. While previous literature postulated QVH, QAH, 
LAF and QSH to be all supporting number of edge states 
growing with the number of layers N (i.e. Z-type in the 
language of classification), we found that the symmetry 
of LAF only protects odd number of edge modes for each 
valley. Hence LAF is a Z2-topological insulator much 
like quantum spin Hall insulator [35]. Transition between 
these different SPT's require closing and re-opening of 
the bulks gap as already been observed in Ref. [34]. 

We predict the naturally occurring tilt boundary [6-8] 
to be the first topological structural defect hosting topo- 



logically protected gapless mode of transport. Most im- 
portantly, our findings on tilt boundaries combined with 
the recent observations [6-8] solve the long standing mys- 
tery of sub-gap transport [1]. Our explanation can be 
tested through proposed transport, scanning tunneling 
spectroscopy and thermopower imaging experiments, and 
scanning SQUID experiments. Experimental confirma- 
tion of the tilt-boundary transport origin of the sub-gap 
transport will open doors to control the sub-gap trans- 
port and enable device application of gated multi-layer 
graphene systems. 

Note added. After completion of this work, a comple- 
mentary preprint [45] , which covers material closely re- 
lated material, has appeared. 
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